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Liquid  metals  in  the  small  radial  gap«  between  the  rot-KS  and  stators  of  homopolar  machines 
represent  low- resistance  electric  current  collectors.  De  i(n  predictions  for  these  liquid-metal 
current  collectors  require  a  thorough  knowledge  of  Hqcd-mctal  flow's  in  a  narrow  gap  between 
a  fixed  and  a  moving  surface,  with  a  strong  applied  magnetic  field  r.nd  a  free  surface  beyond 
each  end  of  the  gap.  The  radial  and  axial  velocities  in  the  secondary  flow  are  reduced  by  a 
strong  axial  or  radial  magnetic  field.  For  a  sufficiently  strong  «kJa.  <hc  azimuthal  momentum 
transport  hy  the  secondary  flow  can  be  neglected.  This  assurer  ion  .etiu'es  the  problem  for  ihc 
primary  azimuthal  velocity  to  a  fully  developed  magnctohydrodynamic  duct  flow  prcbkm 
with  a  moving  wall  and  two  free  surfaces.  Asymptotic  solutions  for  lsi£*  SUrtmann  numbers 
arc  presented  for  skewed  magnetic  fields  with  both  radul  and  axial  compooci'.ss.  Collects 
without  any  electrical  insulation  or  with  insulation  on  the  stator  sides,  or  rotor  sides,  or  both 
are  considered.  Solutions  for  a  putcly  axial  magnetic  field  and  arbitrary  Hartmann  number* 
are  also  presented. 
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I.  INTRODUCTION 

In  a  homopolar  generator,  the  rotation  of  electrically 
conducting  disks  in  a  strong,  axial  magnetic  field  produces  a 
large  dc  electric  current.  This  current  can  be  used  to  drive  a 
homopolar  motor  which  has  a  different  applied  axial  mag¬ 
netic  field  strength,  a  different  number  of  disks,  and  a  differ¬ 
ent  disk  radius.  This  moior-generator  combination  provides 
a  highly  efficient  and  extremely  flexible  method  of  speed 
reduction,  along  with  other  advantages  such  as  maximum 
torque  at  zero  motor  speed.  If  homopolat  devices  can  be 
perfected,  they  wdl  have  a  number  of  important  applica¬ 
tions.1 

The  dc  electric  current,  with  a  typical  current  density  of 
10’  A /m  ’.  must  flow  between  the  lip  of  each  rotating  disk 
(rotor)  and  a  static  current  collector  (stator)  with  a  mini¬ 
mum  voltage  drop.  State  of  the  art  solid  brush-slip  ring  sys¬ 
tems/  such  us  silver  graphite  brushes  on  a  copper  slip  ring 
with  a  protective  humidified  carbon  dioxide  atmosphere, 
typically  have  contact  resistances  in  the  range  of  200-600 
/ill  at  current  densities  un  to  6-7  •  10**  A/m\  For  larger 
current  densities,  the  Joulean  heating  at  the  interface  causes 
excessive  wear.  Furthermore,  because  of  the  required  cool¬ 
ing  and  mechanical  loadings  on  the  brush,  it  is  difficult  to 
achieve  packing  densities  of  the  brush  on  the  sl:p  ring  of 
greater  than  50%-60%.  On  the  other  hand,  liquid  metals,  in 
addition  to  providing  uniform  slip  ring  coverage  and  essen¬ 
tially  wear-free  operation,  have  contact  resistances  of  1.4 /til 
at  7  x  10'*  A/m2.  Furthermore,  liquid-metal  current  collec¬ 
tors  have  been  experimentally  demonstrated  to  run  for  ex¬ 
tended  periods  at  an  order  of  magnitude  larger  current  den¬ 
sity/ 

Liquid-metal  current  collector  design  predictions  re¬ 
quire  a  thorough  understanding  of  the  liquid-metal  flow  in  a 
narrow  gap  between  a  moving  perfect  conductor  and  a  static 
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perfect  conductot  in  the  presence  of  a  strong  magnetic  field 
and  with  a  fret  surface  beyond  each  end  of  the  gap.  One 
objective  is  to  predict  the  voltage  difference  between  the  sta¬ 
tor  and  rotor  for  a  given  electric  current  between  them.  A 
second  objective  is  to  predict  and  minimize  the  viscous  dissi¬ 
pation  and  Joulean  power  losses  in  the  liquid  metal.  A  third 
objective  is  to  insure  that  the  liquid  metal  stays  in  the  radial 
gap  at  all  operating  speeds.  In  present  prototypes,  an  insta¬ 
bility  occurr  at  some  critical  rotation  rate  and  leads  to  ejec¬ 
tion  of  the  liquid  metal  from  the  gap.  This  ejection  appears  to 
involve  a  Kelvin-Helmholtz  instability  at  the  free  surface 
due  to  the  difference  between  the  primary  azimuthal  veloc¬ 
ities  in  the  liquid  metal  and  cover  gas.4  The  axial  gaps  be¬ 
tween  the  sides  of  the  rotor  and  confining  radial  walls  are 
filled  with  an  inert  cover  gas  from  the  shaft  to  the  liquid- 
metal  free  surface  near  the  rotor  tip.  Accurate  predictions  of 
azimuthal  velocities  in  the  liquid  metal  under  various  oper¬ 
ating  conditions  are  needed  to  develop  stable  designs. 

A  typical  value  of  the  magnetic  Reynolds  number 
R„  ^  firoUL  is  0.027.  Here,  fir  *  4jtX  10  *  H/m  is  the 
magnetic  permeability,  a  =  2.38  x  10*’  S/m  is  the  electrical 
conductivity  of  a  liquid  sodium-potassium  eutectic  mixture 
at  31 1  K,  (/a  OR  =s  90  m/sec  is  the  rotor  tip  velocity, 
L  ss  10" 4  m  is  the  radial  gap  between  the  rotor  tip  and  the 
stator,  11  =  300  rad/sec  is  the  angular  velocity,  and  R  a  0.3 
m  is  the  rotor  radius.  Therefore,  we  neglect  the  magnetic 
field  produced  by  the  electric  currents  in  the  liquid  metal. 
The  magnetic  field  in  the  liquid-metal  region  consists  of  the 
radial  and  axial  magnetic  field  components,  which  are  pro¬ 
duced  by  superconducting  coils  around  the  shaft,  and  the 
azimuthal  magnetic  field  component,  which  is  produced  by 
the  electric  currents  in  the  solid  conductors,  Tlte  magnetic 
field  varies  spatially  over  distances  which  are  comparable  to 
R.  Since  L  =  3.3  X 10-4A  and  a  typical  axial  dimension  of 
the  liquid-metal  region  is  26JL  =  8.7  X 10"  3A,  we  neglect  the 
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the  electromagnetic  body  force  opposing  the  secoxtary  flow 
gives  a  characteristk  velocity  for  I  he  radial  and  axial  veloc¬ 
ities,  namely 


FIO  I  L»*<MHl'mei»lrvKtHt«M'*«nlh<slaliH4mJnH>K  IhctoiHtlttMlc 
Md  >««fihv  »(«  mxnuli/ot  by  the  imImI  **p  hn^ih  / 


variations  of  the  magnetic  Held  components  in  the  small  liq- 
uid-metai  region. 

The  assumption  that  t.<B  also  means  that  variations  of 
the  radial  coordinate  arc  negligible,  so  that  we  can  use  the 
Cartesian  coordinate  system  shown  in  Fig  I .  The  origin  is  at 
the  center  of  the  stator,  the  ;  axis  is  parallel  to  the  shaft 
center  line,  the  y  axis  points  radially  inward  toward  thecen- 
ter  line,  and  the  x  axis  is  in  the  a/imuthal  direction.  The 
coordinates  and  lengths  arc  normalized  by  the  radial  gap 
length  /..  so  that  2b/.  is  the  axial  length  of  the  rotor  tip,  and 
each  axial  gap  has  a  length  of  oL  We  use  the  typical  values 
b  m  10  and  a  -  3.  so  that  the  axial  dimension  of  the  liquid- 
metal  region  is  2b/..  Here,  -  If,  and  B  represent  the  radial 
and  axial  components  of  the  magnetic  Ikld,  respectively, 
while  6  is  the  angle  between  (heir  resultant  and  they  axis  and 
B„ «  (fl;  +  /?•)''•'  Wc  need  only  consider  ihc  range 
0«;  tfs.  90*  because  switching  ihc  sign  of  B,  simply  reflects  the 
solution  about  the  z  ^  0  plane,  while  switching  the  signs  of 
both  Br  and  B ,  simply  changes  the  signs  of  certain  variables. 

If  the  shaft  is  horizontal,  then  gravity  produces  a  devi¬ 
ation  from  an  uxisymmclric  frce-surfacc  geometry  and  miIii- 
tion*.  A  simple  application  of  lubrication  theory  indicates 
that  a  dimensionless  parameter  that  reflects  the  relative  devi- 
ation  from  an  axisynimctric  solution  is 

G  e  /#fl ’/.  S/p  U, 

where  /i  K67  kg/m  ‘and  /i  6.5  •  10  ‘  kg/in  sec  are  the 
density  and  viscosity  of  the  liquid  metal,  while#  4,  XI  m/ 
sec’.  A  typical  value  for  (!  is  0.013.  Therefore,  we  neglect 
gravitational  effects  and  treat  only  axisynmielric  flows  in 
which  all  variables  are  independent  of  x.  While  gravity  pro¬ 
duces  only  a  small  perturbation  of  Ihc  axisymnietric  flow, 
this  small  perturbation  may  play  a  key  role  in  determining 
the  ejection  point  in  the  Kclvin-Hclmholtz  instability.  In  a 
future  paper,  we  will  present  a  perturbation  solution  for 
gravitational  effects  for  G'<  I. 

While  the  primary  azimuthal  velocity  u  is  always  com¬ 
parable  to  the  rotor  tip  velocity  U  =  HR,  a  strong  magnetic 
Held  suppresses  the  secondary  flow  involving  the  radial  and 
axial  velocities  —  u  and  w,  respectively.'  A  balance  between 
the  radial  centrifugal  force  due  to  the  azimuthal  velocity  and 


U.  ^ptVB/oB*. 

The  ratio  of  the  characteristic  electromagnetic  body  force 
that  acts  on  the  primary  (lew,  otlRi},,  to  the  characterittk 
transport  of  azimuthal  momentum  by  the  secondary  Dow, 


is  the  interaction  parameter 
<V  -  trB'L/pHVR 


pUJXB 


with  a  typical  value  of  36.2  for  4,  m  6  T.  Here,  we  assume 
that  N>  I,  so  that  the  transport  of  azimuthal  momentum  hy 
the  secondary  flow  is  neglected.  With  this  assumption,  the 
primary  flow  is  decoupled  from  the  secondary  flow  and  we 
can  determine  a  independent  of  a  and  w.  Once  u  k  known,  o 
and  to  can  he  determined,  where  the  centrifugal  body  force 
-  /w’A  'fr.ia  a  known  driving  force  for  the  secondary  flow 
and  fr  ha  unit  vector.  Here,  we  only  present  solutions  for  the 
primary  flow. 

The  assumption  A>  I  is  probably  reasonably  accurate 
for  B„  -  6  T.  but  would  not  be  appropriate  for  machines 
with  weaker  magnetk  (kids.  In  a  future  paper,  we  will  pres¬ 
ent  numerical  results  for  arbitrary  N,  with  full  coupling 
between  the  primary  and  secondary  flows.  The  present 
large-, V  analytical  solutions  and  the  future  artatrary-.V  nu¬ 
merical  solutions  are  complementary.  The  analytical  solu¬ 
tions  arc  sufficiently  simple  to  treat  a  wide  variety  of  operat¬ 
ing  conditions  and  gain  physical  insights  into  the  flow 
phenomena.  The  numerical  solutions,  using  the  approach  of 
Langlois  and  Walker.*1  are  more  accurate,  particularly  for 
weaker  magnetk  fields,  hut  require  sufficient  computer  linn 
such  that  only  n  few  important  cases  can  be  treated  thor¬ 
oughly. 

Since  here  wc  only  consider  the  primary  flow,  the  azi¬ 
muthal  magnetic  field  B,  plays  no  role  in  the  solution.  The 
electric  currents  do  interact  with  B,  to  produce  pressure 
variations  whkh  influence  the  frce-surfacc  locations,  hut 
I  hese  effects  arc  part  of  t  he  secondary  flow  probkm.  Here  wc 
lake  the  frce-surfacc  locations  as  known. 

The  Hartmann  number 


M  BjAn/p)11' 

has  a  typical  value  of  36.3  for  B„  6  T.  The  characteristic 
ratio  of  the  electromagnetic  body  force,  acting  on  ihc  pri¬ 
mary  flow,  1 .1  the  viscous  force  is  M \  For  A/>  I,  the  liquid- 
metal  region  can  be  divided  into  the  subregions  shown  in 
Fig.  2  for  axial  free  surfaces  aligned  with  the  rotor  tip,  i.e..  for 
/,(x)  -  =  1.  Here  G-C3  are  the  :nviscid  core  regions, 

HI-HK  are  the  Hartmann  layers  with  0(M  ')  thickness 
whkh  separate  the  core  regions  from  solid  or  free  surfaces, 
and  FI  and  F2  are  the  interior  or  free  shear  layers  with 
0(M  l/:)  thickness  which  lie  along  the  iwo  magnetic  field 
lines  through  the  comers  where  the  perfectly  conducting  ro¬ 
tor  tip  meets  the  electrically  insulating  free  surface.  Asymp- 
totk  treatments  of  fully  developed  liquid-metal  magnetohy¬ 
drodynamic  (MHD)  duct  flows  for  I  have  been 
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ducting  stator  normalized  by  UBJ*.  Equation  (la)  is  the 
azimuthal  component  of  the  momentum  equation  for  axi- 
symmetric  (low,  with  negligible  azimuthal  momentum 
transport  by  the  secondary  how;  Eqa.  ( lb)  and  ( Ic)  are  the 
if  and  f  components  of  Ohm's  law;  and  Eq.  ( Id)  is  the  con¬ 
servation  of  current  equation. 

Since  d  it  the  dimensionless  deviation  from  the  electric 
potential  in  the  stator,  the  boundary  conditions  at  all  three 
stator  surfaces  are 


I  KS  2  $**fe<ryvo«>  tif  lit*  mtul  for  ,W>  I  mm!  tS*  C<u\oum  «*wvt«. 

iMto  r  *.  v.i  1.  with  «h«  V  [muI kl  10  the  m«t<wOe  Mi  W  lit*  *  *  cw*M 
pUnc 


prevented  in  a  number  of  papers  published  during  the  IWs 
and  were  thoroughly  reviewed  by  Hunt  and  Sherchf.’  The 
prevent  paper  applies  (his  well-established  approach  to  a 
novel  problem  of  technological  importance. 

The  solution  for  the  geometry  shown  in  Fig.  2  for  M >  I 
and  (K0%aretan(«)  -  72’  is  presented  in  Sec.  II.  This  solu- 
lion  for  these  special  fret-surface  locations  illustrates  all  the 
important  elements  of  a  large-, W  asymptotic  solution.  In  Sec. 
III.  wc  show  that  the  Urge-.W  asymptotic  analysis  is  not  ap¬ 
propriate  for  nearly  axial  magnetic  fields  and  wc  present  a 
separate.  arbtlrary-.W  solution  for  the  axial  magnetic  Add 
case  0  90*.  In  See.  IV,  wc  extend  the  analysis  to  other  free- 
surface  locations./,  I  amJ/.  >  I,  with.W>  I  and0<9O*.  In 
Sec.  V,  wc  extend  the  analysts  further  to  include  the  cases 
where  the  sides  of  the  rotor  at  2  ±b  or  the  sid  .'the 
stator  at  c  -  t-  (a  *■  b)  arc  coveted  with  thin  U.  r  an 
electrical  insulator.  Since St  ,  J  *  0.17  for  B„ **  v  .snot 
particularly  small  and  becomes  larger  for  weaker  Helds,  the 
large-, W  solutions  presented  here  are  only  valid  for  very 
strong  magnetic  fields,  A  future  paper  will  present  numeri¬ 
cal  solutions  for  arbitrary  A’ and  St. 

II.  AXIAL  FREE  SURFACES  AT  THE  ROTOR  TIP  WITH  A 
SKEWED  MAGNETIC  FIELD 

Here  wc  prevent  the  asymptotic  solution  for  M >  I  for 
0 *  0-  arctan(tf)  and/,  k/j  I.  It  is  convenient  to  use  a 
different  set  of  Cartesian  coordinates  (x,i/t£)  which  have 
been  rotated  by  an  angle  ((about  the x  axis,  so  that  the  if  axis 
is  parallel  to  the  resultant  of  the  radial  and  axial  magnetic 
field  components,  as  shown  in  Fig.  2.  The  dimensionless  gov¬ 
erning  equations  for  the  primary  flow  are 


./«/'«  ,i‘u\ 

(1*) 

J"  ihf  ’ 

(lb) 

(1c) 

dl,  d$ 

(Id) 

where  u  is  the  azimuthal  velocity  component  normalized  by 
U,j„  and  y.  are  two  components  of  the  electric  current  den¬ 
sity  normalized  by  aUBw  and  d  is  the  deviation  of  the  elec¬ 
tric  potential  function  from  its  value  in  the  perfectly  con- 


u  &  0,  (2a) 

d  v  0.  (2b) 

Since  the  rotor  is  a  perfect  conductor  moving  through  a  mag¬ 
netic  (kid,  the  boundary  conditions  at  the  rotor  face  are 

tnl,  (3a) 

d *4i  + mii 0  »?>♦(£).  (3b) 

where  dn  U  the  dimensionless  potential  at  the  center  of  the 
rotor  face,  i.e.,  at  y  ^  I.  z  0,  or  if  ^  cos  9.  £  —  -  sin  0. 
The  dimensional  voltage  difference  between  the  rotor  and 
stator  at  their  midplanc  *  «  0.  is  d„  UBJL.  We  assume  that 
the  cover  gas  is  an  electrical  insuUtor  with  negligible  shear 
stresses,  so  that  the  boundary  conditions  at  both  free  sur¬ 
faces  arc 


0. 


(4a) 


(3c) 


4-0.  (4h) 

In  the  three  core  regions  CI-C3,  all  derivatives  are 

0(  1),  so  that  the  solution  of  Eqs.  (I)  is 

Ju  ~  0-  <5») 

d,  «d.  -tf*.  (3b) 

j* 

neglecting  0(M  ’)  terms.  Here,  d,  (f)  and /*.(£)  are  inte¬ 
gration  functions  which  are  determined  by  matching  the 
Hartmann  layer  solutions.  The  Hartmann  layers  adjacent  to 
the  stator  or  rotor  surfaces  satisfy  the  boundary  conditions 
(2a)  and  (3a)  and  match  any  core  velocity.  The  Hartmann 
layer  structure  that  accommodates  a  jump  in  the  0(  I )  u 
leads  to  an  0( St  1 )  jump  in  the  electric  potential  across  the 
layer  for  any  perfectly  conducting  wall  whose  surface  is  not 
perpendicular  or  parallel  to  the  magnetic  field  vector.  Since 
the  maximum  jump  in  d  is  0(  St  ’),  theconditions  (2b)  and 
(3b)  can  be  applied  directly  to  the  0(  I )  core  solutions.  For 
the  Hartmann  layers  adjacent  to  the  two  free  surfaces,  the 
0(  1 )  w  is  continuous  through  the  layer.  The  boundary  layer 
structure  develops  if  the  0(  1 )  value  of  du/dy  in  the  core  is 
not  zero  at  y  =  1.  This  structure  involves  a  jump  in  the 
0{M  ~ ')  velocity  and  the  0{M  -I)  normal  electric  current 
density  jy.  Therefore,  the  boundary  condition  (4b)  can  be 
applied  to  the  0(1)  and  0(M~l)  variables  in  the  core  re¬ 
gions  Cl  and  C3. 

For  the  core  regions  Cl  and  C3,  the  solution  (3)  that 
satisfies  the  Hartmann  layer  jump  conditions  is 


j(c  —  0> 

;*•  =  o, 


(6a) 

(6b) 
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4.-0,  (6c) 

y,  “0,  (6d) 

neglecting  0(3/  °)  terms.  The  boundary  condition  (2b)  » 
applied  to  the  0(1)  and  0(3/  1 )  potential*  4  because  there 
is  no  jump  in  0(1)  v  across  any  of  the  surer  Hartmann 
layers  H4,  HJ,  H7,  and  Ht.  In  the  core  C2,  the  0(  I )  van* 
able*  are 


Jri  «  -cosW.  (7a) 

4.i  -  (*cc*ff-£Mnff)’l'.  (7b) 

at,  *  ijccsff- sin0(*  HI  £),  (7c) 

where  the  rotor-surface  potential  ♦(£)  is  defined  by  Eq. 
(3b).  The  0(3/  '  *)  variates  m  core  C2  are  easily  obtained 
by  matching  the  jump  condition*  for  the  Hanmann  layers 
H2  and  H6,  but  these  variables  arc  not  presented  here  be* 
cause  they  represent  minor  corrections,  the  core  solution 
(7)  holds  in  the  region  defined  by 

0<iy  —  £un  disced,  for  )f  i-  sin 0 J  <•  h eos 6. 


The  core  regions  Cl  and  C3  arc  stagnant,  current-free 
regions.  The  core  C2  involves  an  0(  I )  current  density  which 
is  parallel  to  the  magnetic  field  and  equal  to  the  local  poten¬ 
tial  difference  between  the  stator  and  rotor  ♦(  £ )  U'VKtcd  by 
the  distance  between  the  stator  and  rotor  along  a  magnetic 
Add  line  sec  9.  The  Arst  term  in  the  velocity  solution  (7c) 
represents  linear  Couctic  flow  ami  the  second  term  is  an 
additional  velocity  associated  with  the  current  flow.  This 
additional  velocity  is  independent  of »/,  linear  in  f.  decreases 
lincady  as  d„  is  ‘ncrcascd.  and  vanishes  for  a  purely  radial 
magnetic  Add  6 0.  Wc  discuss  the  core  solutions  further 
after  we  present  the  free  shear  layer  solutions. 

The  core  solutions  (6c)  and  (7b)  indicate  that  each  free 
shear  layer  must  accommodate  a  jump  in  the  0( I )  potential 
4-  The  magnitude  of  this  jump  increases  linearly  from  7.cro  at 
the  tutor  to  4\  m  4u  -  b  cos  0  of  4i  «  4u  +  A  cos  0  at  the 
rotor  comer  for  the  free  shear  layer  FI  or  F2,  respectively. 
With  an  0(  I )  potential  4>  Eq.  ( lb)  implies  that  JH  is  0(  I )  as 
wdl.  Since  d/d$  »  0(3/  nJ)  inside  a  free  shear  layer,  Eq. 
( Id)  implies  that  jl  u  0(3/  m).  Equation  ( Ic)  then  im¬ 
plies  that  v  is  very  large,  namely  0(3/  us) .  The  Arst  and  last 
terms  in  Eq.  ( la)  now  balance,  which  is  why  the  free  shear 
layer  thickness  is  0(3/  ,,1)>  Wc  consider  only  FI  because 
iu  solution  reduces  to  a  universal  solution  which  applies  for 
layer  F2  and  any  similar  free  shear  layer. 

For  the  free  shear  layer  FI,  wc  rescale  the  dependent 


and  independent  variables  by  substituting 

4°*4i4,U,£)  +  0{M  "*).  (8a) 

u  =  3/,,V,A  ,,lMl-f)  +  0(l).  (*b) 

y,  *  4*1  ’  %■(!,$)  +  0(3/  •  •«),  (8c) 

it  =3/- m4,X  J%(/,£)  +  0(3/  - '),  (Id) 
f=  (7  +  3sin0 +Un0sin0),l  \  (8e) 

f  =  3/  ,w(£  -f  sin  0  +  6  cos  0)X  ~ ,n,  (80 


where  X  =  sec  6  is  the  length  ofthe  free  shear  layer  from  the 
sUtor  to  the  rotor  comer.  The  coordinate  /  is  a  local  rescal¬ 
ing  of  the  17  coordinate,  so  that  /  =  0  at  the  sUtor  and  /  =  1  at 
the  rotor  comer.  The  coordinate  £  is  the  stretched  local  £ 


coordinate,  so  that  the  free  shear  layer  solution  matches  the 
Cl  or  C2  core  solution  as  £-  -  «e  or  £-  «,  respectively. 
With  the  factors  involving  4,  and  X  in  Eqs.  (I),  the  free 
shear  layer  variables  4r>  and  jtK  are  universal  func¬ 

tions  of  the  rescaled  coordinates  (/,£)  and  are  completely 
independent  of  the  parameters  4*0, «,  and  b.  With  the  sub¬ 
stitutions  (t),  the  leading  terms  m  Eqs.  ( I )  become 


<9a) 


(%) 


(9c) 


(9d) 


B4> 

J*  ’  ”  *  * 

94r 

*»#  «  -rr  - 

*J4> 

ft*  ’ 

There  are  Hartmann  layers  with  0(3/  ')  thickness 
which  separate  the  free  shear  layer  from  the  stator,  rotor, 
and  free  surface.  The  jump  in  4  across  the  I  lartmann  layers 
adjacent  to  the  stator  and  rotor  is.* ’gain  0(3/  '),sothat  the 
boundary  conditions  (2b)  and  (3b)  can  be  applied  to  the 
0(  I )  fret  shear  layer  4 ■  The  jump  in  jf  across  the  frec-sur- 
face  Hartmann  layer  is  0(3/  *  \  so  that  the  boundary  con¬ 
dition  (4b)  can  be  applied  to  the  0(  I )  free  shear  layer  cur¬ 
rent  density.  These  conditions  become 

4t  -  0,  at  /  -  0,  for  -«<*£<«, 

4,  -  I.  at  l  -■  I,  for  Oe£eao. 


(lOu) 

(10b) 


^L~0. 

dr 


(10c) 


Mia) 

(lib) 


at  r*  1,  for  -  00  <£<Q. 

Matching  the  core  solutions  (6c)  and  (7b)  gives 
4,  -*0,  as  £-•  -  00, 

4)  -7,  as  £  —  ae . 

Once  wc  And  the  solution  4r  that  satisfies  Eq.  (9d)  and  the 
conditions  ( 10)  and  ( 1 1 ).  the  other  variables  are  given  by 
Eqs.  (9a)-(9c). 

The  separation  of  variables  solution  of  Eq.  (9d)  for 
-  *  *  £-  0.  which  satisfies  conditions  (10a),  (10c),  and 
(I  la),  is 

X  Nin(2«i,/)c‘**'  (.-f„  cos(o„£) 


4, 


for 


a 


(12) 


t  Bm  sin(«-.£)|. 

wlierc 

a„  |J(2m  +  l)tr|,'\ 

The  corresponding  separation  of  variables  solution  for 
0«;£<w,  which  satisfies  conditions  (lOu),  (10b),  and 
(lib),  is 


4t-i+  X  **"1**0*  "'‘(C,  cos(fl,£) 


(13) 


+  Dn  sin(/?.£)  J,  for  0<£<  oc, 

where 

0m  =  (Jn»r)'». 

The  constants  Am,  Bm,  C„  and  Dm  are  determined  by 
the  four  conditions  that  4,-  and  its  Arst  three  derivatives  with 
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respect  to  £  are  continuous  at£»0.  Hie  governing  equation 
(9d)  then  ensures  that  all  derivatives  with  respect  io  £  are 
continuous  at  £  --  0.  Since  the  Fourier  scries  in  the  solutions 
C 12)  and  ( 13)  are  not  orthogonal  to  each  other,  we  truncate 
both  scries  and  define  a  residual  from  the  conditions  at 
f*0: 


(i,o  )  -  uo*)j  di. 


where  the  truncated  solutions  ( 12)  and  (13)  arc  used  to 
compute  values  at  £  0  and  at  £  -  0  * .  respectively.  This 
residual  is  minimized  with  rc'pcet  to  the  constants  .4„, 

C„,  and  /),  and  this  minimization  gives  a  set  of  simulta- 
neous,  linear,  algebraic  equations  For  these  constants.  These 
equations  arc  solved  with  Gauss  elimination  and  the  truncat¬ 
ed  senes  ( 12)  and  (13)  arc  used  to  compute  the  other  Free 
shear  layer  variables  From  Fqs  (9a)-(9c). 

The  profiles  of  u,  at  r  -  0.2, 0.4, 0.6. 0.8.  0.94,  and  1.0 
arc  presented  in  Fig.  3.  At  r»l.0,  (i)  ut  *0  For 
-  6.  £<  2  9,  (ii)  ut  increases  rapidly  From  zero  at 

£.*  -  2.9to0.7at£  -  0.  and  (iii)  ut  has  a  discontinuity  at 
£■  0  aFter  which  u,  ■-  0for£t*Q.  For£<  -0.1,  the  pro* 
files  at  r »  0.94  and  1.0  are  very  close.  However,  the  discon* 
tmuity  at  £  -  0  For  /  -  1.0  is  replaced  by  a  rapid  decrease 
between  f  ‘  -  -  0.1  and  £  -  0.9  For  f  -  0.94.  At  /  v  0.94,  u,. 
«  0  For  0.9.  £.  2.  As  t  decreases  From  0.94,  the  peak  veloc¬ 
ity  decreases  to  0.423  at  £  -  0.44  For  r »  0.8.  to  0.278  at 

£  «  -  0.5  For  r «  0.6.  to  0.174  at  £  -  -  0.53  For  / »  0.4, 
and  to  0.084  at  £  -  0.56  For  i  ~  0.2.  In  addition,  the  ve¬ 

locity  profiles  spread  out  as  r  decreases.  Equation  (9c)  and 
conditions  (II)  show  that 


r. 


mi  that  the  Dow  earned  by  each  oF  the  profiles  in  Fig.  3  is 
proportional  to  f. 

Walker  er  al.*M  treat  a  similar  Free  r  Sear  layer  with  the 
present  Free  surFacc  replaced  by  a  solid  electrical  insulator. 
Their  boundary  value  problem  is  the  same  as  the  present 


FIG.  J  Free  shear  layer  velocity  u,  al  r  k  0.2, 0.4, 0.6. 0.8, 0.94,  and  1.0. 


problem  (9)-(ll)  except  that  the  Tree-surface  condition 
( 10c)  is  replaced  by 

— —  *  ?-■ ,  at  t  m  I,  For  -oo<£«0,  (14) 

m 

which  is  the  I  lartmann  layer  matching  condition  at  the  solid 
insulator.  This  difference  precludes  the  use  r>F  the  present 
separation  oF  variables  solutions.  With  a  solid  insulator,  a 
Fourier  transForm  reduces  the  problem  to  an  integral  equa¬ 
tion  which  must  be  solved  numerically.'*  In  spite  oF  the  Fact 
that  the  problem  with  a  solid  insulator  requires  a  much  more 
complex  solution  technique,  the  results  For  the  velocity  pro¬ 
files  and  the  other  variables  arc  very  similar  to  the  present 
results. 

The  solution  in  the  Free  shear  layer  F2  in  terms  of  the 


same  universal  Functions  u,  ,  4, ,  and;’.,,  is 

4  *  d.d,  (r.£)  +  OlM  ,  J).  (15a) 

v*  -  M '  }dwi  '  •'«, (/.£)  +  0(1).  ( 15b) 

4?>  %  (r.£)  +OIM  "*').  (15c) 

Jfm  -  M  +  OC.V  ').  (I5d) 

r  --  (if  +  sin  0tan0  -  6 sin  (i)A  ',  (I5e) 

£*  -.¥"*(£  hsinfl-heosM  "i.  (I5F) 


Here,  the  potential  jump  across  the  Free  shear  layer  increases 
linearly  from  zero  at  the  stator  to  4S «  +  b  cos  0  althe 

rotor  corner,  while  X  sec  0  is  the  length  of  the  layer  From 
the  stator  to  the  rotor. 

For  no  net  electric  current  between  the  stator  and  rotor, 
^u0,  4t  -s  -6cos0,  and  4.  *6  cos  A  The  solutions 
(8b)  and  ( I  Sb)  indicate  that  the  Free  shear  layers  are  identi¬ 
cal  jets  in  the  -  .e  direction  For  this  case.  The  rotor  tip  and 
liquid  metal  in  the  Free  shear  layers  are  going  in  opposite 
azimuthal  directions.  The  velocity  at  the  free  surface  very 
near  each  rotor  corner  is  approximately 
-  U.7.W 1  ,}6(cos  0)  *'  \  As  4»  is  increased  from  zero  to  drive 
a  net  current  From  the  rotor  to  the  stator,  the  Free  shear  layer 
F2  becomes  a  stronger  jet  in  the  -  .r  direction,  while  layer 
FI  decreases  in  strength  to  zero  ffow  at  4»  *  6  cos  0.  For 
4„  >  b  cos  0,  layer  FI  becomes  a  jet  in  the  +  x  direction.  As 
4, ,  is  decreased  From  zero  to  drive  a  net  current  From  the 
stator  to  rotor,  the  above  description  of  layer  FI  applies  to 
layer  F2  and  vice  versa. 

For  core  C2,  the  net  dimensionless  electric  current  from 
the  rotor  to  the  stator,  per  unit  A.t.  is 

24 J)  cos'  U.  (16) 

Equation  (16)  is  multiplied  by  2xuUB,d.R  to  obtain  the  net 
dimensional  current  From  the  rotor  to  the  stator.  As  0  in¬ 
creases  From  zero,  the  electrical  resistance  of  the  liquid-met¬ 
al  region  increases  as  sec*  0:  This  is  because  (he  length  of  the 
current  lines  increases  as  sec  0,  while  the  area  perpendicular 
to  the  current  direction  decreases  as  cos  0.  The  core  regions 
C I  and  C3  carry  no  current.  The  two  Free  shear  layers  allow 
a  small  additional  current  to  ffow  between  the  stator  and 
rotor  because  electric  current  lines  are  allowed  to  Fringe  an 
0(<V  l,:)  distance  beyond  the  lines  £=  —  sin0±6cos0. 
However,  the  ratio  of  this  additional  Free  shear  layer  current 
totheC2corecurrentiscomparableto4f  ~  >l2b  ~  '(sec  0)in. 
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TW*  ratio  s*  quite  small,  except  as  0  approaches  9 0*.  which  is 
the  ease  treated  in  See.  III.  While  the  current-line  fringing 
inside  the  free  shear  layers  docs  not  significantly  increase  the 
net  current,  it  docs  produce  a  small  0{M  "*)  electric  cur* 
rent  component  perpendicular  to  the  magnetic  field.  The 
electromagnetic  body  force  doc  to  this  small  jt  drives  the  jets 
with  0(M  us)  u  inside  the  free  shear  layers. 

Here  we  present  t  he  <?{  At  * )  and  Of  .W '  * )  dimensionless 
power  losses  due  to  Jog  lean  heating  and  viscous  dissipation, 
per  unit  A.e  Total  dimensional  power  losses  arc  obtained  hy 
multiplyingcach  result  by  2xR/tV*.  ThconlyOf.W  *)  dimen¬ 
sionless  power  loss  is  the  Joutcan  heating  in  the  core  C2  This 
loss  is  given  by  the  integral  of  the  square  of  solution  (7a) 
over  the  core  region  and  the  result  is 

l.W*'cos</Cd,‘  d,‘}  (17) 

The  OiM '  '“)  contributions  to  the  power  losses  arc  (i)  the 
viscous  dissipation  due  to  the  large  gradient  of  the  large  vc* 
loc,.y  in  each  free  shear  layer  ami  (it)  (he  Joukan  healing 
due  to  the  fringing  of  th  ♦  0(  I )  dee  me  current  density  insnk 
the  free  shear  layers.  With  the  solutions  (8)  and  ( 15),  these 
power  losses  arc  given  hy  coefficients  times  integrals  of 
(ihif  /dS)*  and  j ^  over  the  region  0-  /-  I.  a;  •  <T  *  * 
Since  these  universal  functions  arc  independent  of  all  param* 
eters,  these  integrals  give  sums  of  the  known  coefficients  in 
the  separation  of  variables  solutions  { 12)  and  ( H)  These 
sums  give  specific  numerical  values  which  apply  to  all  simt* 
lar  free  shear  layers.  The  results  feu  the  dimensionless  sis* 
cons  dissipation  and  Joukan  heating  ar«~ 

OmV'tfM  '  (IXa) 

0.407.W  ’  tyA  '  \  flkh) 

respectively,  where  d,  d,  or  d  for  layer  I- 1  or  1-2 

III.  AXIAL  MAGNETIC  FIELD 

The  solutions  in  Sec  II  apply  forU-  U  arclantu)  7’" 
fore  “  3.  Fora  purely  axial  magnetic  field,  0  ‘XT,  H,  t). 
and  R,  R„  and  wc  use  the  original  (.e.r.z)  Cartesian  co* 
ordinates  shown  hi  Fig  2  With  the  same  iioodimcnMonali* 
ration,  Eqs,  ( I J  arc  replaced  hy 


0  j,  1  •'(  ‘(*2  ♦  ,? 

"  \,h  •  1:  J 

(19a) 

j  ,  ** 

J>  tl;  * 

(1%) 

,/d 

''  *  “• 

( 19c) 

*L+*LW0. 

0y  Oz 

(I9d) 

The  (low  is  now  symmetric  about  ihez 

0  plane,  so  (hat  we 

need  only  solve  Eqs.  ( 19)  for  0-  y ■  1  and  0-  :■  u  +  A.  The 

boundary  conditions  arc 

u  =  0,  at  y  -  0, 

(20a) 

d  *  0.  »t  y  «  0, 

(20b,' 

~  =  0,  at  z  ~  0, 
oz 

(20c) 

£  =  0,  at  z  =  0, 
oz 

(20d) 

kbO,  at  z m a  +  A, 

(20e) 

d  -  0,  at  z  *  •  +  A, 

(200 

u  ^l,  »t/w|,  for  Ode  A. 

(20g) 

d  d,H  at  /  *  1,  fcr0s*<A. 

(20*3 

•=  0,  at  y  as  |,  for  b<z<k*  +  A, 

Of 

(20») 

j,  -  0,  at  jr  -  1,  for  A*  x«*u  f  A. 

(20$) 

Since  the  magnetic  field  is  now  parallel  to  the  rotor  surface, 
this  surface  is  at  a  smgk  dimensionless  elect  he  potential  dw 

In  the  asymptotic  solution  for  ,V>  I,  u  -  d  ** Jf  **  J, 
0  to  all  orders  m  ,V,  except  in  a  thin  boundary  Uycr  which 
hasan0(*W  ,  ?)  thickness  and  is  adjacent  to  the  free  and 
rotor  surfaces  at y  -  I.  All  the  electric  current  from  the  sta¬ 
tor  to  the  rotor  (lows  axially  inside  this  Uycr  near  the  free 
surfaces  and  enters  the  rotor  near  the  corners.  The  azimuthal 
velocity  m  is  OlM '  ’}  inside  this  Uycr.  In  other  words,  as 
0  -90\  the  two  free  shear  layers  from  the  rotor  corners 
merge  to  become  a  frcc-xurfacc,  rotor-face  boundary  Uycr, 
the  current-carrying  core  C2  disappears;  and  the  stagnant, 
current-free  core  C3  expands  to  fill  most  of  the  liquid-metal 
region  However,  the  Urge-.W  asymptotic  solution  is  not  ap¬ 
propriate  for  an  axial  or  nearly  axial  magnetic  field  in  homo- 
polar  devices. 

Figure  J  shows  that  the  jc»s  in  the  free  shear  Uycrs  arc 
confined  to  \  I,  so  that  Eq,  (8f)  indicates  that  the 
free  shear  layer  thickness  IS'U  1  ?.W  1 '  f'or  II  arctan(u), 
the  length  of  each  free  shear  layer  is  A  -  sec  0  The  ratio  of 
the  free  shear  Uycr  thickness  to  the  -  2Aeos0ofeorcC2 
is 

2Csec0)*  M  1  'b  1 

The  Urge-.W  asymptotic  analysis  is  approprUte  for 
U  U  arctan  («)  as  long  as  the  above  rat  hi  is  small:  For 
ll  -15".  V  16  3.  and  h  10.  the  ratio  »s  0  056.  F'or  live 
axial  field  casc.X  becomes  the  entire  knglh  of  the  0\M  *'*) 
tlrvk  ness  layer  adjacent  to  y  I,  mi  that  A  2ta  i  b )  As¬ 
suming  AC  is  still  4.  the  raim  of  the  A>*  of  lhi>  layer  to  the 
entire  Ay  I  is 

4|2(«  *  6)/,Vp  \ 

this  ratio  must  be  small  for  the  Urgc-.W  asymptotic  solution 
to  be  valid  for  the  axial  field  case:  For  a  3.  b  -  10,  and 
.W  36.3  this  ratio  is  3.39.  Therefore,  viscous  efforts  are 
significant  everywhere  and  wc  must  l real  ,\f  as  an  0{  I )  pa¬ 
rameter 

Wc  introduce  ilr-  electric  current  sirvumfonclion 
lily.:),  whK'h  is  defined  by 

J.  M  1  *£.  (21a) 

<fc 

j,  M  (21b) 

Ry 

Sincey,  and  j,  arc  even  and  odd  functions  ofz,  respectively,  A 
is  an  odd  function  of z,  We  eliminate  d  from  Eqs.  ( 19b)  and 
( 19c),  so  that  Eqs.  ( 19)  become 


+  M—~  0, 
* 


(22*) 
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djr  dr  <fc 

TV  boundary  conditions  (20a),  (20c),  (20e),  (20*),  and 
(20i)  on  a  still  apply,  but  the  conditions  (20b),  (20d), 


(200*  (20h),  and  (2Qj)  on  4  *re replaced  by 

~  «  0,  at  /  ■  0. 

to 

A  -  0,  at  x  0, 

(23a) 

(23b) 

—  * 0.  at  tstf  +  4. 

<« 

(23c) 

“»0,  at  y*  1,  for  (Kx<A. 

dr 

A**  A*  atywl,  for  A«xsa  +  A. 

(23d) 

(23e) 

The  constant  A„  is  proportional  to  the  total  dimensional  elec¬ 
tric  current  from  the  stator  to  the  rotor,  which  is 


4irA,,V  ,tf(/MX.»4aAa(//((^v)',^ 

Once  the  solutions  for  u  snd  A  are  known,  we  introduce  them 
into  Eq.  (19c)  and  integrate  with  respect  to  /  to  obtain 
dOgr)  and  +» 

Equations  (22)  and  the  boundary  conditions  (20a), 
(20c),  (20c), and  (23a)-(23c),exceptaty »  I,  are  satisfied 
by  the  separation  of  variables  solutions 

a  -  -  £  costf.x)  [<?,A„  4.  M„  (fu  , 

(24a) 

A  «  j^sintf.x)  [c.f,.  +  //.  -  “-^)] .  (24b) 


where  <7,  and  //.  are  arbitrary  constants  and 

(2n+  l)r/2(tf  +  A),  (24c) 

m,-(2)’,,,«,(0  +  A/,d.  ,),/,+  l)”a.  (24d) 

m,  « (2)  w,d«((l  +  A/Sd,  *),/*  —  (24e) 

A',,  «  stnh(m,y)sin(m,y),  (240 

A’,.  -  cosh(m,y)sin(m)y),  (24g) 

A*,,  -  sinh(m,/)cos(m,y),  (24h) 

A«.  -  cosh(m,y)co»(m,y).  (24i) 


The  constants  G „  and  Hm  are  determined  by  the  remaining 
boundary  conditions  (20g),  (20i),  (23d),  and  (23e)  at 
y  •*  I.  We  truncate  the  scries  and  define  a  residual 

-ri-'Msn* 

♦n($j,+,*-w,i*  <25' 

where  the  integrands  are  given  by  the  series  (24)  evaluated 
aty  =s  I.  When  the  residual  (25)  is  minimiaed  with  respect 
to  the  constants  Gm  and  H„,  it  gives  a  set  of  simultaneous, 
linear,  algebraic  equations  for  these  constants.  These  equa¬ 
tions  are  solved  using  Gauss  elimination. 

Some  typical  results  for  a  *  3,  A  —  10,  M  *40,  and 
Ao=  ±  10  are  presented  in  Fig.  4.  They  scale  is  stretched  by 
a  fretor  of  5  relative  to  the  x  scale.  Figures  4(a)  and  4(c) 


FIG,  4,  Elect  nc  cumtM  Ivu*  A  m  e mm  and  velocity  awgaitads  lim* 
«  m  comt  for  m atioi  magnetic Md  *Hli *■  ).  km  10, aod  M ■  40.  (a) 
Cocmrt  Imwv  10.  tb)  velocity  weyiMtttdci  foe  A,  »  10.  (c)>  WM 

iMctforA,*  -  Iftwtd  (d)  velocity  fer  A, "  -  *9. 


show  that  electric  current  prefers  to  Bow  axially  along  the 
magnetic  Held  lines  near  the  free  surface,  rather  than  flow 
across  the  field  lines  in  the  central  part  of  the  radial  gap.  For 
A**  10  or  -  10, 78%  or  60%  of  the  total  current  leaves  or 
enters  the  stator  side  at  x  »  a  +  A,  rather  than  the  axial  sur¬ 
face  at  /  =  0.  For  current  from  the  stator  to  the  rotor,  Fig. 
4(b)  shows  that  the  maximum  velocity  near  the 

rotor  corner.  For  current  from  the  rotor  to  the  stator,  Fig. 
4(d)  shows  that  the  minimum  velocity  at  the  free  surface  is 
lie  —  1.4  near  the  rotor  corner,  but  the  minimum  velocity 
is  v  =  -2nearyst0.5for0<x<6. 

In  a  homopolar  device,  there  arc  two  groups  of  rotors  or 
disks  which  are  separated  axially  along  the  shaft.  Each  disk 
in  one  group  is  electrically  connected  along  the  shaft  with 
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one  disk  in  (he  other  group.  Electric  current  flows  radially 
inward  inside  all  the  disks  in  one  group  and  radially  outward 
inside  all  the  disks  in  the  other  group.  The  magnet  coils  cre¬ 
ate  asial  magnetic  fields  in  opposite  directions  for  the  two 
groups  of  disks.  Therefore,  the  electromagnetic  body  force  is 
in  the  same  azimuthal  direction  in  all  the  disks  in  both 
groups.  The  analysis  here  assumes  that  the  disk's  angular 
velocity  vector  and  the  asial  magnetic  field  *e  ;,n  the  same 
direction,  i.c.,  that  U  ^  11  Al  >  0  and  Bt  >  0.  For  the  present 
solution,  if  we  replace  B,  -  B(t  with  B,  -  B„  and  we  re¬ 
place  A  by  -  A.  then  the  velocities  arc  exactly  the  same.  In 
addition,  the  elect  rie  currents  arc  exactly  the  same,  bat  arc  in 
the  opposite  direction.  The  results  for  A„  >  0  apply  for  cur¬ 
rent  from  the  stator  to  the  rotor  for  B,  ->  0  and  for  current 
from  the  rotor  to  the  stator  for  B,  <0  and  vice  versa  for 
Ad  Therefore,  A„  **  0  corresponds  to  all  liquid-metal  cur¬ 
rent  collectors  in  a  homopolar  motor.  In  a  motor,  the  elec¬ 
tromagnetic  body  force  in  the  disks  and  liquid  metal  is  in  the 
direction  of  rotation,  so  that  u  is  increased  to  u  >  I ,  as  showr. 
in  Fig.  '4(b).  On  the  other  hand,  An>  6  corresponds  to  all 
current  collectors  in  a  homopolar  generator.  In  a  generator, 
the  azimuthal  electromagnetic  body  force  in  the  disks  and 
liquid  metal  is  in  the  opposite  azimuthal  direction  from  the 
disk  rotation,  so  that  u  is  decreased  to  u-0.  as  shown  in  Fig. 
4(d). 

The  value  of  d„  is  a  linear  function  of  A„.  For  example, 
foi  u  -  3,  A  10,  and  M  10. 

A,  -  -  0.456*0.082  43Air. 

The  total  Joulean  healing  and  viscous  dissipation  arc  para¬ 
bolic  functions  of  A,,.  For  example,  for  u  -  3,  b  ■■  10,  and 
M 10.  the  dimensionless  total  Joulean  heating  and  viscous 
dissipation  arc 

2.736  i  0  954A,,  t  Q.2909A‘, 

21.65  0.984A,,  +  1.J43A'. 

respectis  ety.  Total  dimensional  power  losses  arc  obtained  by 
multiplying  each  result  by  2 zK/tU'' 


IV.  GENERAL  FREE-SURE  ACE  LOCATIONS  ANO  A 
STRONG,  SKEWED  MAGNETIC  FIELD 

Here  we  present  the  extension  of  the  large-Af  asymptotic 
analysis  for  a  skewed  magnetic  field  to  general  frce-surfacc 
locations  at  y «/,(x)  for  -  «  A%x*  -  Aandaty  -  /.(z) 
for  At.zso  +  A,  as  shown  in  Fig.  5.  The  solution  (6)  still 
applies  in  core  regions  Cl  and  C3.  The  electric  potential  of 
all  rotor  surfaces  is  ♦(£)  =  dn  +  sin  0  +  £.  In  core  C2,  the 
solution  in  the  radial  gap  is  still  given  by  the  results  (7)  for 
|£  +  sin  0|  <Acos  9.  For  —  Acor^— A  ~  A)sin  0<$ 
<  —  A  cos  B  -  sin  dr  is  a  linear  function  of  »/,  which 
equals  zero  at  the  stator  and  ♦(£)  at  the  rotor,  whiley'^.  and 
uc  are  given  by  Lqs.  ( lb)  and  (5c),  respectively.  There  is  a 
weak  free  shear  layer  at  £  =  —  sin  6  —  A  cos  8  along  the 
magnetic  field  line  through  the  corner  at  y  =  1,  z  =  -  A. 
The  electric  potential  is  continuous  across  this  layer,  but  u , 
is  discontinuous  between  the  two  parts  of  core  C2.  The  free 
shear  layer  structure  that  matches  this  jump  in  ue  involves 
only  an  0(  1 )  velocity,  so  that  it  is  a  much  weaker  layer  than 


I'KJ  5  Sutoqtiw*  I  He  Hqufcl  meut  far  My  I  gcM-ul  fm-Mufaec 
H<*Ihv»s  7 He  lUttttunn  b)*»\  *n»J  a  wwl  flee  'Hear  Layer  al 

l  Mit  II  ictKitOxnon 


layers  FI  and  F2,  which  involve  an  velocity.  In 

core  C4, 

d.  ♦tf).  «.  I.  U  - h .  -0.  (26) 

neglecting  0(.U  **)  terms. 

The  solution  in  the  free  shear  layer  FI  is  given  by  expres¬ 
sions  (8) terms  of  the  same  universal  functions  d, .  ut . 
and  y,# ,  but  now  X  is  the  free  shear  layer  length  from  the 
stator  to  the  point  where  the  free  surface  and  rotor  side  meet 
aty  /,v  h),:*  -  A.  while 

d,  d(l  r  sin 0 1 1  -/,{  A))  Aco$0 

is  at  tha  point. 

For  the  free  shear  layer  F2,  we  again  rcscak  the  depen¬ 
dent  and  independent  variables  by  introducing  expressions 
(15)  with  X  -  see  0  and  dj  ‘  d.t  +  A  cos  0.  However,  the 
functions  d, ,  u,  ,  and  ju  are  not  the  same  functions  as 
before;  thus  we  add  the  subscript  2  to  denote  these  new  func¬ 
tions.  These  rescaled  variables  are  functions  of  (/.£)  for 
0  -  «  cfe  so,  where  e  equals  the  length  of  the  free 

shear  layer  F2  from  the  stator  to  the  free  surface  divided  by 
see  0.  The  functions  dM,  ,,  and  jin  depend  on  one 

parameter  «r,  so  that  they  are  no  longer  universal  functions. 
There  is  a  slot  in  the  domain  at  t  --  I  for  ()«  £<•  <*  respesent- 
ing  the  two  surfaces  ot  'he  rotor,  so  that  the  variables  need 
not  be  continuous  at  t  -  I  for  £V0.  Equations  (9)  and  the 
boundary  conditions  (10a),  (IQb),  and  (I la)  still  apply. 
The  condition  ( 10b)  applies  for  both  sides  of  the  slot  at 
l  I .  The  Tree-surface  condition  ( 10c )  becomes 

«7d 

-  0,  at  r  =  ar,  for  -«<£<<»,  (27) 

at 

while  matching  core  regions  C2  and  C4  gives 

d/i-/.  as  f-oo,  forOsKl,  (28a) 

d«-l.  as  f—  oo,  for  K/<*r.  (28b) 

The  separation  of  variables  solution  (12)  applies  for 
-  o©  <£v,0  and  0  «./<;*■,  with  am  replaced  by 

nm  =  l  (2/m  +  I)jt/4*]i,j. 

The  separation  of  variables  solution  (13)  applies  for 
<  oo  and  0<r<  1.  For  0<f  <  oo  and  we  intro¬ 

duce 
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*„«!  +  i»in[2 riu-'))'  n{ 

*  -a 

X  (/\  co*(y»£)  +  0,  Mn(n?)].  (29) 

where 

>'»  -  |(2 k  »  I )  s/4  (  a'  -  U|,#* 

The  coefficients  A„,  /?„,  /*t,  and  &  arc  determined 

by  minimizing  the  tame  residual  E  used  in  Sec.  II,  except 
that  the  limns  *.)f  integration  are  from  r  -  0  to  r  -  a\ 

Specific  results  for  this  new  one-parameter  free  shear 
layerarc  not  presented  here.  ForCk/s,  I,  «hc  velocity  profile* 
arc  very  similar  to  those  tn  Fig.  J.  As «  increases  from  I,  the 
velocity  profiles  spread  out  and  the  peak  velocity  decreases, 
but  with  'he  constraint  that 
* 

ii,  ,  I.  for  all  I 

a* 

As  the  left  free  surface  rises  fromy  —  I  loy -*/i(z),the 
free  shear  layer  FI  simply  follows  the  point  y  */,(  -  b), 
z  -  -  b  and  the  free-surface  velocity  is  still  zero,  except 
near  the  rotor  side.  The  only  changes  in  the  free  shear  layer 
FI  arises  from  changes  in<i  and  d,.  Some  additional  current 
flows  between  the  stator  and  rotor  through  the  new  pari  of 
core  C2.  As  the  right  free  surface  rises  from  y  **  1  to 
y  v /,(z),  a  new  core  C4  with  u  m  I  develops  and  ihc  free 
shear  layer  F2  intersects  the  free  surface.  The  free-surface 
velocity  is  I  from  the  rotor  side  to  the  free  shear  layer, 
±  0(M ,  :)  inside  the  free  shear  layer,  and  zero  from  the 
free  shear  layer  to  the  stator  side.  The  volume  flux  per  unit 
Ai?  in  the  free  shear  layer  F2  increases  lineatly  from  zero  at 
the  stator  to  a  value  at  the  rotor  comer  and  then  keeps  that 
constant  value  from  the  rotor  corner  to  the  free  surface. 


V.  INSULATING  LAYERS  ON  THE  SIDES  OF  THE  ROTOR 
OR  STATOR 

Here  we  present  the  extension  of  the  analysis  of  Sec.  IV 
todevices  with  thin  layers  of  electrical  insulators  on  the  sides 
of  the  rotor  at  :  s-  t  b  or  on  the  side  of  the  stator  at 
z  -  i :  (<»  +  b).  The  thickness  of  the  electrical  insulators  is 
assumed  to  be  much  smaller  than  M  ’  Ui,  so  that  they  are 
thinner  than  the  free  shear  layers. 

The  subregions  of  the  liquid  metnl  for  insulators  on  the 
sides  of  the  rotor  arc  shown  in  Fig.  6.  The  solution  (6), 
neglecting  0[M  j  terms,  again  holds  in  core  regions  Cl 
and  C.V  The  Hartmann  layer  udjacent  to  the  insulator  at 
z  b  must  match  the  jump  in  u  from  zero  in  core  Cl  to  I 
at  the  rotor  surface.  This  Hartmann  layer  involves  an  0(  1 ) 
current  density  j,  and  the  circuit  for  this  current  must  be 
completed  through  a  free  shear  layer  at 
£=•  -  b  cos  0  -  /,  (  —  A)sin  0,  i.e„  along  the  magnetic  field 
line  through  the  point  where  the  free  surface  and  insulator 
met*'  This  implies  an  0(M  u2)  current  density  jH  inside 
this  tree  shear  layer,  which  in  turn  Implies  an  0(1)  fluid 
velocity  u.  Thus  this  free  shear  layer  is  another  weak  layer 
involving  only  an  0(  1 )  velocity. 

The  core  region  C2  occupies  the  same  region  and  has  the 
same  solution  (7),  neglecting  0(M  “')  terms,  as  in  Sec.  II. 
In  the  core  C4, 
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l-’lti  6  SubfcgHVis  uf  Ihc  Itqukl  n<ct»l  for  My  t  »<*d  thin  ck<lfic»l  mmiU- 
tors  im  ihc  Mdo  of  the  inter  at  i  -  *  A  The  Hartman*  layers  and  »  weak 
fret  shoir  layer  at  f  -  At(H0/,f  -  AJun  d  arc  shown 


ut  ■*  I,  d,  *s  C  -t-  £  -  b  coa  0  +  sin  $, 

neglecting  OlM  !)  terms.  The  constant  Cis  determined  by 
matching  the  free  shear  layer  F2. 

The  rescaling  (8)  and  the  boundary  value  problem  (9)- 
(II)  for  the  free  shear  layer  FI  is  the  same  as  that  presented 
in  Sec.  II,  except  that  the  free-surface  condition  ( ICAc)  is 
replaced  by  the  solid  insulator  condition  (14).  Walker  tl 
aiy  present  the  solution  for  this  modified  free  shear  layer: 
Although  the  analysis  is  much  more  complex,  the  results  arc 
surprisingly  similar  to  those  presented  in  Fig.  3. 

The  boundary  value  problem  for  the  free  shear  layer  F2 
is  the  same  as  that  presented  in  Sec.  IV,  except  that  the 
boundary  condition  on  the  top  of  the  slot  in  the  domain  be¬ 
comes 

+  at  f-r,  for  0<f < «  (30) 

and  matching  core  C4  gives 

df:~C,  as  £~oo.  for  l<r<*.  (31) 

The  change  from  t‘te  boundary  condition  ( 10b)  at  /  *  I  *  to 
condition  (30)  precludes  the  relatively  simple  separation  of 
variables  solution  discussed  in  Sec.  IV.  Without  solving  for 
this  new  layer,  we  can  draw  some  conclusions  about  its  solu¬ 
tion.  The  integral  of  un  from  -  oo  tof  =  »  equals  the 

jump  in  >t  each  value  of  f.  For  0</<  I,  this  integral 
equals  f,  sc  that  the  velocity  profiles  here  cannot  deviate 
significantly  '")m  those  in  Fig.  3.  With  an  insulator  at  z  =  b, 
there  is  no  ct.  c-*;j  potential  jump  for  1  <K*,  but  the  jump 
here  must  be  independent  of  t,  i.e.,  it  must  equal  C.  The  jet  for 
1  <  i<,k  and  the  associated  potential  jump  C  are  driven  by 
continuity  of  the  electric  potential  and  the  electric  cur¬ 
rent  density —  —  dfF/dt  at  /=  1  for  -  oo  <f  <0.  The 
velocity  profiles  may  spread  out  as  /  increases  from  I ,  but  the 
total  flow  per  unit  A /  in  each  profile  must  be  the  same. 

The  effect  of  adding  an  insulating  layer  on  the  side  of  the 
rotor  at  z  =  —  b  is  to  move  the  free  shear  layer  FI  from  the 
magnetic  field  line  through  the  point y  =/,(  —  A),z  =  —b 
to  the  field  line  through  the  pointy  =  1,  z  =  —b.  The  strip 
of  fluid  between  bese  two  lines  becomes  stagnant  and  cur¬ 
rent  free,  except  in  the  Hartmann  layer  at  z  =  —  b.  The  net 
electric  current  between  the  stator  and  rotor  is  once  again 
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2 4xfi  cot5  9,  as  in  Sec.  II.  The  only  effect  of  adding  ihe  insu- 
lating  layer  at  z  ■  b  ii  t  slightly  alter  the  velocity  profiles 
inside  the  free  shear  layer  F2. 

Increasing  the  magnitude  of/,(z)  in  Fig.  6  does  not 
change  anything.  However,  if/,(:)  is  increased  sufficiently 
so  that  the  free  shear  layer  F2  intersects  the  stator  rather 
than  the  Tree  surface,  then  several  changes  occur  in  core  04 
and  the  free  shear  layer  F2  ThcC4corcsolutionwithu,  **  I 
still  applies  in  the  triangle  formed  by  the  insulator  at :  -  b, 
thtf  free  surface  at  y  «/.(:),  and  the  magnetic  field  line 
through  the  point  where  the  free  surface  and  stator  meet  at 
y  */jU  +  b),  :  a  »  b .  There  is  a  new  core  05  between 
this  magnetic  field  lincand  the  free  shear  layer  F2,  i.c.,  in  the 
strip 

(«  +  b) cos  0  -/j(n  +  A)sin  0<-£\6  cos  0  -  sin  0 . 

(52) 

In  this  new  core,  ut  -  4,  -  j*  j(,  -  0.  neglecting 
0{\t  :)  terms.  There  is  a  weak  free  shear  layer  between 
cores  04  and  05.  This  weak  layer  matches  the  jump  in  u , 
from  zero  in  05  to  I  in  04  and  completes  the  circuit  for  0(  I ) 
current  density  j,  in  the  Hartmann  layer  between  the  insula* 
tor  at  z  *  b  and  core  05.  In  the  (x>undary  value  problem  for 
the  free  shear  layer  F2,  conditions  (27)  and  (51)  are  re* 
placed  by 

4n  “0.  at  /  -  a",  for  -  «  » •£>  «,  (55a) 

4n  asf~«,  for  I  <  /nat.  (55b) 

Again,  the  velocity  profiles  for Q  j<  I  must  closely  resemble 
those  presented  in  Fig,  5,  but  now  the  integral  of  u,  from 
f  »  -  «  to  f  »  no  must  equal  zero  for  1  *;  /*,a\  The  profile 
at  /  a  1  *  consists  of  -  u,  for  f  >  0  and  l  u,  for  f -■  0.  by 
continuity  with  the  solution  at  /  I  .  As  t  increases  from  I, 
this  basic  character  persists,  but  the  magnitude  of  the  veloc¬ 
ity  quickly  decreases  to  zero.  Therefore,  for  /  •  I,  there  is  a 
double  jet  with  no  net  flow  which  is  concentrated  very  near 
the  rotor  corner. 

We  now  consider  the  addition  of  insulating  layers  on  (he 
sides  of  the  stator  at :  ~  y  la  h  b)  to  the  situations  shown 
in  Figs.  5  and  6,  i.c.,  with  insulators  only  on  the  stator  sides 
and  on  both  the  stator  and  rotor  sides,  respectively.  For  Figs. 
5  and  6  wc  have  assumed  that  /,(  b )  is  sufficiently  small 
such  that  the  magnetic  field  line  through  the  point 
y  */,(  -  b),  :  -  -  b  intersects  the  stator  at  y  0  and 
file  +  h )  is  sufficiently  small  such  that  the  free  shear  inter¬ 
sects  the  free  surface.  For  these  low  frce-surface  Ideations, 
the  addition  of  insulators  at  r  t  ( a  i  h)  has  absolutely 
no  effect  on  (he  solutions  with  or  without  insulators  at 
z—±b. 

lf/,(  -  b)  is  increased,  so  that  the  magnetic  field  line 
through  the  pointy -/,(  -b),z~  -  b  intersects  the  insu¬ 
lated  stator  side  at  z  -  -  a  -  b,  then  a  new  core  C6  occurs 
between  this  magnetic  field  lincand  the  field  line  through  the 
corner y  =  0,  z  =  —a  -  b,  i.c.,  in  the  strip 

—  b cos  0  —  /,(  -  6)sin  0r£»,  -  (a  +  6)cos  0.  (34) 

The  Cl  solution,  with  u,  =  0,  still  applies  in  the  triangle 
above  £  =  —  b  cos  0 — /,  (  -  fc)sin  0. 

If  the  side  of  the  rotor  at  z  =  -  b  is  not  insulated,  then 
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u,  m  i,  4,  *  vi$),  «  o,  (35) 

neglecting  0(M  *})  terms,  in  the  new  core  C6.  The  free 
shear  layer  FI  lies  along  the  magnetic  field  line  through  the 
comer  at  y  ■  0,  z  »  -  a  -  b  and  the  C2  core  solution  ap¬ 
plies  to  the  right  of  the  free  shear  layer.  The  solution  of 
Walker  *r  at.**  applies  to  layer  FI.  The  principal  change  is 
that  the  potential  jump  across  this  layer  is  now  zero  at  the 
rotor  surface  at  t  ®  I  and^,  *  do  +  sin  9  -  (a  +  b) cm  0at 
the  stator  comer  at  /  »  0.  The  insulating  free  surface  has 
moved  above  layer  FI  and  this  layer  now  stems  from  Ihe 
junction  of  the  insulating  and  perfectly  conducting  stator 
walls.  Therefore,  the  velocity  profiles  closely  resemble  those 
in  Fig.  3  with  /replaced  by  (I  -  /).  Moving  in  the  +y  direc¬ 
tion  from  the  free  shear  layer  FI,  u(  »  i  in  the  new  core  C6; 
u  decreases  from  1  to  zero  inside  the  weak  free  shear  layer 
along  the  magnetic  field  line  through  the  pointy  j/,(  -  b), 
z  -  b‘,  and  a,  -  0  in  core  C5  to  the  free  surface. 

We  now  consider  the  problem  lor  a  high  left  free  surface 
with  insulators  at  both  :  ah  and  :  b.  The  free 

shear  layer  FI  is  m  the  position  shown  »n  Fig.  6  and  its  solu¬ 
tion  is  unchanged  by  the  addition  of  the  insulator  at 
:  -  ~  a  -  A.  The  Cl  core  solution  with  u,  -«  4,  y„ 

-  0.  neglecting  OlM  })  terms,  applies  in  thestrip  between 
the  free  shear  layer  FI  and  the  magnetic  field  line  through 
the  stator  corner  at  y  -  0,  z  -  a  ~  b.  The  new  core  C6 
occupies  the  strip  defined  by  the  limits  (54)  and  has  a  solu¬ 
tion 

«.  0,5,  4.  -  0.5|£  »•  («  + A)cos0),  (56) 

withyv  y4,  0,  neglecting  OlM  3 )  terms.  By  conserva¬ 

tion  of  current,  the/,,  inside  the  Hartmann  layers  adjacent  to 
the  insulators' at:  -  -  6  and  at:  -  u  -  b  must  be  equal 
and  opposite.  Since  the  jump  in  u  across  these  Hartmann 
layers  is  proportional  to  these  currents,  u,  is  0.5  between  a 
moving  and  a  fixed  insulator.  In  ihclriangularcorcCI  above 
the  magnetic  field  line  through  the  point  y  */,(  -  b ), 
b. 

4,  *  0.5a  cos  0  -  0.5/,  (  -  6)sin  0, 

with  a,  */,„  l,  0,  neglecting  OlM  5)  terms.  Moving 
upward  from  layer  FI,  u  is  first  zero;  then  0.5;  and  finally, 
zero  again. 

With  a  high  right  free  surface,  the  triangular  core  C4 
extends  to  the  stator  at :  ~  a  y  b  and  there  is  a  new  core  CS 
between  core  C4  and  the  free  layer  1**2,  i.c.,  in  the  strip  de¬ 
fined  by  the  limits  (32).  In  the  cores  C3  and  C4,  u,  ~  0  and 
I,  respectively,  with  or  without  insulators  on  either  side.  If 
the  stator  side  at  z  -  u  +  b  is  insulated,  but  the  rotor  side  at 
:  -  b  is  not  insulated,  then  the  solution  (35)  applies  in  the 
new  core  C5.  Therefore,  u,  ~  I  throughout  the  axial  gap 
above  the  free  shear  layer  F2.  The  only  distinction  between 
cores  C4  and  C5  is  that  there  is  a  weak  free  shear  layer  be¬ 
tween  them  which  completes  the  circuit  for  the  jt  in  the 
Hartmann  layer  between  the  new  core  C5  and  the  insulator 
at  z  =  a  +  b.  The  boundary  value  problem  for  the  free  shear 
layer  F2  is  the  same  as  that  discussed  in  Sec.  IV,  except  that 
the  boundary  condition  (27)  is  replaced  by 

34n  024fi  .  .  r  f  tm\ 

— — =■■■■■■.-,  at  /  =  *■,  for  -  oo  <f<  oo.(37) 
at  ag 
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This  change  precludes  the  separation  of  variables  solutions 
presented  in  Sec.  IV,  but  it  does  not  change  the  fact  that  the 
integral  of  ur  from  -  »  tof  a  «e  is  I  for  all  1<K*. 

If  both  sides  at  x  »  6  and  x  »  fl  +  b  are  insulated,  then 
thesolution(36)  with^,  =  0.5£  +  C  +  sin  cos  0  ap¬ 
plies  in  the  new  core  C5.  The  constant  C  is  determined  by 
matching  the  free  shear  layer  F2.  As  we  move  up  the  axial 
gap  from  layer  F2,  u,  first  equals  0.S  in  core  C3  and  then 
equals  1  in  the  core  C4.  The  boundary  value  problem  for  the 
free  shear  layer  F2  is  the  same  as  that  discussed  earlier  in  this 
section,  except  that  the  boundary  condition  (33a)  is  re¬ 
placed  by  condition  (37)  and  the  matching  condition  (31 ) 
applies.  Again,  the  integral  of  i from  fa*  —  «  to  f  >■  co 
equals  a  constant  C  for  all  1  </<*,  where  C  is  different  from 
that  for  an  insulated  rotor  side  and  a  perfectly  conducting 
stator  side  and  is  determined  by  solving  the  free  shear  layer 
problem. 


VI.  CONCLUSIONS 

This  paper  presents  solutions  for  the  primary  azimuthal 
flow  in  liquid-metal  current  collectors  for  homopoUr  de¬ 
vices.  If  the  local  magnetic  field  is  strong  and  has  both  radial 
and  axial  components,  then  an  asymptotic  analysis  for  large 
Hartmann  numbers  is  appropriate.  The  solutions  in  the  in- 
viscid  core  regions  are  very  simple.  The  solutions  for  the  free 
shear  layers  along  magnetic  field  lines  through  junctions  of 
perfect  conductors  and  insulators  are  either  presented  here 
or  discussed.  This  asymptotic  analysis  leads  to  simple  ex¬ 
pressions  for  the  net  electric  current  and  the  power  losses 
due  to  Joulean  heating  and  viscous  dissipation. 

For  a  nearly  axial  magnetic  field,  the  asymptotic  analy¬ 
sis  is  not  appropriate.  Arbitrary  Hartmann  number  solu¬ 
tions  are  presented  for  an  axial  magnetic  field  and  axial  free 
surfaces  aligned  with  the  rotcr  face  at  /  m  |.  This  analysis 
can  be  extended  to  other  axial  free  surfaces  at  y  »/, 
b/i  =  const.  One  separation  of  variables  solution  is  used 
for(ky<  l,0<x<fi  and  satisfies  the  boundary  conditions 
aty  -  O.y  **  I,  and  z  =  0.A  second  separation  of  variables 
solution  applies  for  0 A<x«r  -f  b  and  satisfies  the 
boundary  conditions  at  y  =  0,y  *s  /3,  and  x  *  a  +  b.  The  re¬ 
sidual  is  formed  from  the  boundary  conditions  at  x  «  b  for 
1  ~A  «nd  the  continuity  of  u,  du/dt,  A,  and  dh/dt  at 
z  ~  b  for  0<y<  I .  This  method  works  for  problems  without 
any  insulating  layers  or  with  insulating  layers  on  the  stator 
sides,  or  rotor  sides,  or  both.  This  method  can  be  extended  to 
different  free-surface  locations,  where /,  and  fj  are  different 
constants.  The  symmetry  about  x  =  0  is  then  lost  and  the 
number  of  coefficients  in  the  separation  of  variables  solu¬ 
tions  is  doubled. 


A  typical  value  of  the  Reynolds  number  for  the  primary 
flow,  Re« pUL/ft,  is  12000.  At  present,  it  is  uncertain 
whether  this  flow  is  laminar  or  turbulent.  An  ordinary 
Couette  flow  with  this  Reynolds  number  would  probably  be 
turbulent.  However,  a  strong  magnetic  field  suppresses  tur¬ 
bulence  and  increases  the  value  of  the  critical  Reynolds 
number.  This  magnetic  suppression  of  turbulence  depends 
on  the  electrical  characteristics  of  the  walls  and  is  strongest 
with  highly  conducting  walls  which  allow  current  circula¬ 
tions  which  damp  the  turbulent  eddies.  In  addition,  the  gap 
is  very  small,  so  that  wall  damping  is  strong.  These  consider¬ 
ations  would  indicate  that  the  flow  is  laminar.  On  the  other 
hand,  MHD  effects  lead  to  liquid-metal  velocities  which 
may  be  as  large  as  3  (/or  may  be  negative,  so  that  Re  based  on 
(/and  L  may  underestimate  the  largest  velocity  gradient  in 
the  flow.  The  velocity  profile  depends  on  the  amount  of  elec¬ 
tric  current  between  the  stator  and  rotor,  so  that  the  flow 
may  be  laminar  for  one  current  and  turbulent  for  another. 
The  laminarization  of  a  turbulent  flow  by  a  strong  magnetic 
field  has  been  studied  experimentally  for  flows  in  circular 
pipes  and  rectangular  ducts,  but  these  results  may  not  apply 
for  the  present  Couette  flows  in  a  very  small  gap  with  strong 
electromagnetic  pumping;  there  is  clearly  a  need  for  appro¬ 
priate  experiments. 
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